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Abstract 

Searches for pair-produced scalar leptoquarks are performed using 20 fb"^ 
of proton-proton collision data provided by the LHC and recorded by the AT¬ 
LAS detector at = 8 TeV. Events with two electrons (muons) and two or 
more jets in the final state are used to search for first (second)-generation 
leptoquarks. The results from two previously published ATLAS analyses are 
interpreted in terms of third-generation leptoquarks decaying to bvjbvr and 
tvjvj final states. No statistically significant excess above the Standard Model 
expectation is observed in any channel and scalar leptoquarks are excluded at 
95% CL with masses up to mLQi < 1050 GeV for first-generation leptoquarks, 
mi^qi < 1000 GeV for second-generation leptoquarks, mLQ 3 < 625 GeV for 
third-generation leptoquarks in the bvrbvr channel, and 200 < mLQ 3 < 640 GeV 
in the channel. 


© 2016 CERN for the benefit of the ATLAS Collaboration. 

Reproduction of this article or parts of it is allowed as specified in the CC-BY-3.0 license. 




Contents 


1 

Introduction 

3 

2 

The ATLAS detector 

4 

3 

Data and Monte Carlo samples 

5 


3.1 

Monte Carlo for background predictions 

5 

4 

Searches for first- and second-generation LQs 

6 


4.1 

Trigger and data collection 

7 


4.2 

Object selection 

7 


4.3 

Event pre-selection 

8 


4.4 

Signal regions 

9 


4.5 

Background estimation 

10 



4.5.1 Control regions for Z/y*-f-jets and tt backgrounds 

10 



4.5.2 Kinematic distributions 

11 


4.6 

Systematic uncertainties 

11 


4.7 

Results 

13 

5 

Search for third-generation LQs in the bv^bv-r channel 

15 


5.1 

Object and event selection 

15 


5.2 

Background estimation 

17 


5.3 

Results 

17 

6 

Search for third-generation LQs in the tVjtVr channel 

19 


6.1 

Object and event selection 

21 


6.2 

Background estimation 

22 


6.3 

Results 

22 

7 

Summary and conclusions 

26 


2 


1 Introduction 


Leptoquarks (LQ) are predicted by many extensions of the Standard Model (SM) [1-7] and may 
provide an explanation for the many observed similarities between the quark and lepton sectors of 
the SM. LQs are colour-triplet bosons with fractional electric charge. They carry non-zero values 
of both baryon and lepton number [8]. They can be either scalar or vector bosons and are expected 
to decay directly to lepton-quark pairs (where the lepton can be either charged or neutral). 

The coupling strength between scalar LQs and the lepton-quark pairs depends on a single Yukawa 
coupling, termed The additional magnetic moment and electric quadrupole moment in¬ 

teractions of vector LQs are governed by two coupling constants [9]. The coupling constants for 
both the scalar and vector LQs, and the branching fraction of the LQ decay into a quark and a 
charged lepton, /?, are model dependent. The production cross-section and couplings of vector LQs 
are enhanced relative to the contribution of scalar LQs, although the acceptance is expected to be 
similar in both cases. This analysis considers the simpler scenario of scalar LQ pair-production, 
for which the form of the interaction is known and which provides more conservative limits on LQ 
pair-production than for vector LQ pair-production. 

In proton-proton collisions, LQs can be produced singly and in pairs. The production of single 
LQs, which happens at hadron colliders in association with a lepton, depends directly on the un¬ 
known Yukawa coupling However, LQ pair-production is not sensitive to the value of the 

coupling. In pp collisions with a centre-of-mass energy = 8 TeV, the dominant pair-production 
mechanism for LQ masses below ~1 TeV is gluon fusion, while the ^^-annihilation production 
process becomes increasingly important with increasing LQ mass. 

The minimal Buchmiiller-Ruckl-Wyler model (mBRW) [10] is used as a benchmark model for 
scalar LQ production. It postulates additional constraints on the LQ properties, namely that the 
couplings have to be purely chiral, and makes the assumption that LQs are grouped into three 
families (first, second and third-generation) that couple only to leptons and quarks within the same 
generation. The latter requirement excludes the possibility of flavour-changing neutral currents 
(FCNC) [11], which have not been observed to date. 

Previous searches for pair-produced LQs have been performed by the ATLAS Collaboration with 
1.03 fb“' of data collected ai ^[s -1 TeV, excluding at 95% confidence level (CL) the existence 
of scalar LQs with masses below 660 (607) GeV for first-generation LQs at = 1 (0.5) [12] and 
685 (594) GeV for second-generation LQs at /3 = 1 (0.5) [13]. The CMS Collaboration excluded 
at 95% CL the existence of scalar LQs with masses below 830 (640) GeV for first-generation LQs 
at = 1 (0.5) and 840 (650) GeV for second-generation LQs at yS = 1 (0.5) with 5.0 fb“^ of data 
collected at = 1 TeV [14]. 

Pair-produced third-generation scalar LQs decaying to bvrbvr have been excluded by the CMS 
Collaboration for masses below 700 GeV at P = 0, and for masses below 560 GeV over the full p 
range using 19.7 fb“^ of data collected at ^[s = % TeV [15]. Third-generation scalar LQs have 
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been excluded in the bT^bT~ channel at jS = 1 for masses up to 740 GeV by the CMS Collaboration 
using 19.7 fb“' of data collected at = 8 TeV [16], and by the ATLAS Collaboration at jS = 1 for 
masses up to 534 GeV using 4.7 fb“' of data collected at ^Js -1 TeV [17]. The CMS Collaboration 
also excluded third-generation scalar LQs in the tT~tT'^ channel at/3 = 1 for masses up to 685 GeV 
using 19.7 fb“^ of data collected at = 8 TeV [15]. 

In this paper, searches for pair-produced first- and second-generation scalar LQs (LQl and LQ2, 
respectively) are performed by selecting events with two electrons or muons plus two jets in the 
final state (denoted by eejj and nidjj, respectively). In addition, limits are placed on pair-produced 
third-generation scalar LQs (LQ3) by reinterpreting ATLAS searches for supersymmetry (SUSY) 
in two different channels [18,19]. LQ production and decay mechanisms can be similar to those 
of stop quarks (?) and sbottom quarks (b). For example, ti —> ttx^x^ gives the same event topology 
as LQ3LQ3 —>tv^?Vr in the limit where the neutralino (y-^) is massless. Two ATLAS analyses 
optimised for these SUSY processes are therefore used to set limits on the equivalent LQ decay 
processes: LQ3LQ3 -^bvrbvr andLQ3LQ3 -^tVrivr- 

The results for each LQ3 channel cannot be combined since the parent LQs have different electric 
charges in the two cases {-^e for the LQ3 LQ3 —>bvrbvr channel and for the LQ3 LQ3 —^tVrtVr 
channel, where e is the elementary electric charge). The branching fractions of LQ3 decays to 
bv-r and tv-^ are assumed to be equal to 100% in each case. Although complementary decays of 
a charge -^e (|e) LQ into a tT~tT^ {bT^bT~) final stale are also allowed, kinematic suppression 
factors which favour LQ decays to /^-quarks over t-quarks and the relative strengths of the Yukawa 
couplings would have to be considered. Since these suppression factors are model dependent, limits 
are not provided as a function of /3 for the LQ3 channels. 


2 The ATLAS detector 


The ATLAS experiment [20] is a multi-purpose detector with a forward-backward symmetric cyl¬ 
indrical geometry and nearly An coverage in solid angle. The three major sub-components of 
ATLAS are the tracking detector, the calorimeter and the muon spectrometer. Charged-particle 
tracks and vertices are reconstructed by the inner detector (ID) tracking system, comprising silicon 
pixel and microstrip detectors covering the pseudorapidity ^ range I 77 I < 2.5, and a straw tube tracker 
that covers \r]\ < 2.0. The ID is immersed in a homogeneous 2 T magnetic field provided by a solen¬ 
oid. Elecfron, photon, jet and tau energies are measured with sampling calorimeters. The ATLAS 
calorimeter system covers a pseudorapidity range of I 77 I < 4.9. Within the region I 77 I < 3.2, elec¬ 
tromagnetic calorimetry is provided by barrel and endcap high-granularity lead/liquid argon (LAr) 


* ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the 
detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the (/-axis 
points upward. Cylindrical coordinates (r, 0) are used in the transverse plane, 0 being the azimuthal angle around the 
beam pipe. The pseudorapidity is defined in terms of the polar angle 6 as ;/ = - In tan(0/2). 
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calorimeters, with an additional thin LAr presampler covering \r]\ < 1.8, to correct for energy loss in 
material upstream of the calorimeters. Hadronic calorimetry is provided by a steel/scintillator-tile 
calorimeter, segmented into three barrel structures within \t]\ < 1.7, and two copper/LAr hadronic 
endcap calorimeters. The forward region (3.1 < \q\ < 4.9) is instrumented by a LAr calorimeter 
with copper (electromagnetic) and tungsten (hadronic) absorbers. Surrounding the calorimeters is 
a muon spectrometer (MS) with air-core toroids, a system of precision tracking chambers providing 
coverage over \r]\ < 2.7, and detectors with triggering capabilities over \r]\ < 2.4 to provide precise 
muon identification and momentum measurements. 


3 Data and Monte Carlo samples 

The results presented here are based on proton-proton collision data at a centre-of-mass energy of 
1/5 = 8 TeV, collected by the ATLAS detector at the LHC during 2012. Data samples correspond¬ 
ing to an integrated luminosity of 20.3 fb“^ are used by all channels except for the LQ3LQ3 ^ 
bvrbvr analysis which uses 20.1 fb“^ because of requirements made by the trigger used in the ana¬ 
lysis. 

Simulated signal events of pair-produced scalar LQs decaying to e^e~qq, n^iTqq, tv^tYr, and 
bvrbvr final slates are produced using the Pythia 8.160 [21] event generator with CTEQ6L1 [22] 
parton distribution functions (PDFs). The coupling which determines the LQ lifetime and 

width [23] is set to VO.Ol x Ana, where or is the fine-structure constant. This value gives the LQ 
a full width of less than 100 MeV, which is much smaller than the detector resolution. For LQ 
masses in the ranges considered here (200 GeV < wlq < 1300 GeV, in steps of 50 GeV), the value 
of the coupling used is such that the LQs can be considered to decay promptly. The production 
cross-section of pair-produced LQs is assumed to be independent of the coupling strength. The 
signal process is normalised to the expected next-to-leading-order (NLO) cross-sections for scalar 
LQ pair-production [24]. The signal production cross-section is 23.5 fb for a LQ mass of 600 GeV, 
and 0.40 fb for a 1 TeV LQ and is the same for each generation. 


3.1 Monte Carlo for background predictions 

The Monte Carlo (MC) samples used to estimate the contributions from SM backgrounds to the 
LQl and LQ2 searches are discussed here. Details about the MC models used for estimating 
backgrounds in the LQ3 searches are available in Refs. [18] (for the bvrbvr channel) and [19] (for 
the tVrtVr channel). 

The MC samples used to model the Z/y*-i-jets background with a dilepton invariant mass {mu) less 
than 120 GeV and high-mass Drell-Yan backgrounds {mu > 120 GeV) are generated with SHERPA 
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1.4.1 [25]. The high-mass Drell-Yan samples are generated assuming massive c- and Zj-quarks 
instead of the conventional massless treatment. 

Samples of tt events are produced with POWHEG box [26, 27] interfaced with PYTHIA 6. MC 
samples representing the WW, fkZ, and ZZ diboson decays are generated with HERWIG 6.52 [28] 
and use the AUET2 [29] values for the tunable parameters (the ‘AUET2 MC tune’). Samples of 
single-top-quark events in the Wt and ^-channel are generated with MC@NLO 4.01 [30,31] and 
the AUET2 MC tune, while the t-channel samples are generated with AcerMC 3.8 [32] interfaced 
with PYTHIA 8 and use the AUET2B [33] MC tune. The hadronisation and parton showering of 
the samples produced with MCONLO are done using HERWIG 6.52 coupled to JIMMY 4.31 [34]. 
The VT-i-jets samples are produced with ALPGEN 2.14 interfaced with JIMMY 4.31, also with 
the AUET2 MC tune applied. The choice of PDFs used to produce the MC simulated samples is 
generator dependent: AcerMC, PYTHIA, HERWIG and ALPGEN use CTEQ6L1, while MC@NL0 uses 
CTIO [35]. For all samples, the detector response is modelled [36] using GEANT4 [37], except for 
the Drell-Yan background samples, which use a fast detector simulation where the calorimeter 
response is parameterised. The differences between fast and full simulation in terms of kinematic 
spectra and modelling of relevant objects are evaluated to be negligible. 

The cross-sections of background processes used in the analysis are taken from theoretical predic¬ 
tions. Single-top production cross-sections in the 5-channel [38], t-channel [39], and in associated 
production with a W boson [40], are calculated to NLO-i-NNLL accuracy. IT-i-jets and Z tt 
cross-sections with NNLO accuracy are used [41]. The cross-sections for WW, WZ, and ZZ pro¬ 
cesses are calculated at NLO [42,43]. The theoretical cross-section for WW production is scaled 
by a factor 1.2 and the uncertainty is increased by an extra 20%, in order to take into account the 
ATLAS [44] and CMS measurements [45], which showed an excess in data at the level of 20% (see 
Refs. [46,47] for more discussion about possible causes of the excess). 

For the Z/y*-i-jets and tt backgrounds, LO and NLO cross-sections, respectively are used. These 
backgrounds are constrained using two control regions (CRs), as described in Section 4.5. 


4 Searches for first- and second-generation LQs 

The first- and second-generation analyses exploit similarities in the final states and use common 
search strategies to select dilepton plus dijet final states. Control regions are used to constrain 
estimates of the dominant backgrounds to the data. A set of discriminating variables is used to 
define signal regions (SRs) that are used for a counting analysis. 
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4.1 Trigger and data collection 


Selected data events are required to have all relevant components of the ATLAS detector in good 
working condition. For the LQl {eejj) analysis, the trigger requires at least two electromagnetic 
calorimeter clusters, defined as energy deposits in the cells of the electromagnetic calorimeter. The 
leading cluster is required to have transverse momentum px > 35 GeV and the sub-leading one 
Pt > 25 GeV. This trigger selects electrons without imposing any requirement on the isolation and 
this allows a data-driven estimate of the background contribution from jets in the final state that 
pass the electron selection, as described in detail in Ref. [48]. The trigger is 98% efficient with 
respect to the offline selection, which requires px above 40 (30) GeV for the leading (sub-leading) 
electron. 

For the LQ2 (pp/y) analysis, events are selected from data using a trigger which requires the pres¬ 
ence of at least one muon candidate in the event with px above 36 GeV. This trigger is fully efflcient 
relative to the offline selection for muons with px above 40 GeV [49]. 

4.2 Object selection 

Electrons are selected and identified by imposing requirements on the shape of the cluster of en¬ 
ergy deposits in the calorimeter, as well as on the quality of the track, and on the track-to-cluster 
matching. The identification efficiency is on average 85% [50]. Electron candidates must have 
transverse energy Ex > 30 GeV and \r]\ < 2.47. Electron candidates associated with clusters in 
the transition region between the barrel and endcap calorimeters (1.37 < \r]\ < 1.52) are excluded. 
All electrons are required to be reconstructed with cluster-based or combined cluster- and track- 
based algorithms and to satisfy calorimeter quality criteria. Requirements are made on the trans¬ 
verse (|r/ol) and longitudinal (|zol) impact parameters of the electron relative to the primary vertex 
and must satisfy |dol < 1 rnm and |zol < 5 mm. In addition, electrons are required to be isol¬ 
ated by imposing requirements on the measured in the calorimeter within a cone of size 

V(^ -I- (Ad>)2 = 0.2 around the electron cluster excluding the electron cluster energy, and 
corrected to account for leakage (i.e. energy deposited by the electron outside of the cluster) and 
the average number of proton-proton interactions per bunch-crossing. The isolation requirements 
are optimised for high-px electrons following the strategy in Ref. [48]. The leading electron is 
required to have < 0.007 x Ex -i- 5 GeV, and the sub-leading electron is required to have 

eAR< 0.2 < Q Q 22 X Ex -t 6 GeV. 

Muon tracks are reconstructed independently in the ID and the MS. Tracks are required to have 
a minimum number of hits in each system, and must be compatible in terms of geometrical and 
momentum matching. In particular, in order to prevent mis-measurements at high px, muons are 
required to have hits in all three MS stations, as described in Ref. [48]. In order to increase the 
muon identification efficiency, when one muon in the event satisfies the three-stations requirement, 
the criteria for the second muon in the event are relaxed to require hits in only two MS stations. 
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Information from both the ID and MS is used in a combined fit to refine the measurement of the 
momentum of each muon [51]. Muon candidates are required to have pj > 40 GeV, \t]\ < 2.4, Ir/ol < 
0.2 mm and |zol <1.0 mm. Muons must also pass a relative-isolation requirement Ipt < 

0.2, where is the sum of the transverse momenta of all the tracks with pj above I GeV 

(except for the muon track) within a cone of A/? < 0.2 around the muon track, and px is the 
transverse momentum of the muon. 

Jets are reconstructed from clusters of energy deposits detected in the calorimeter using the anti- 
kt algorithm [52] with a radius parameter R - 0.4 [53]. They are calibrated using energy- and 
77 -dependent correction factors derived from simulation and with residual corrections from in-situ 
measurements. The jets used in the analysis must satisfy px > 30 GeV and |? 7 | < 2.8. Jets recon¬ 
structed within a cone of AT? = 0.4 around a selected electron or muon are removed. Additional jet 
quality criteria are also applied to remove fake jets caused by detector effects. A detailed descrip¬ 
tion of the jet energy scale measurement and its systematic uncertainties is given in Ref. [54]. 


4.3 Event pre-selection 

Multiple pp interactions during bunch-crossings (pile-up) can give rise to multiple reconstructed 
vertices in events. The primary vertex of the event, from which the leptons are required to originate, 
is defined as the one with the largest sum of squared transverse momenta of its associated tracks. 
Events are selected if they contain a primary vertex with at least three associated tracks satisfying 
PT,track > 0.4 GcV 

MC events are corrected to better describe the data by applying a per-event weight to match the 
distribution of the average number of primary vertices observed in data. A weighting factor is also 
applied in order to improve the modelling of the vertex position in z. Scale factors are applied to 
account for differences in lepton identification and selection efficiency between data and MC simu¬ 
lation. The scale factors depend on the lepton kinematics and are described in detail in Ref. [51] for 
muons, and in Ref. [55] for electrons. The energy and momentum of the selected physics objects 
are corrected to account for the resolution and scale measured in data, as described in Ref. [51] for 
muons, in Ref. [55] for electrons and in Ref. [54] for jets. 

Events are selected in the eejj channel if they contain exactly two electrons with px > 40 (30) GeV 
for the leading (sub-leading) electron and at least two jets with px > 30 GeV. Eor the ppjj channel, 
events are selected if they contain exactly two muons with px > 40 GeV and opposite-sign charge, 
and at least two jets with px > 30 GeV. No requirements are placed on the charges of the electron 
candidates due to inefficiencies in determining the charge of high-px tracks associated with elec¬ 
trons. These sets of requirements form the basic event ‘pre-selection’ for the analyses, which is 
used to build the control and signal regions discussed in the following sections. 


4.4 Signal regions 


After applying the event pre-selection requirements, a set of signal regions is defined using addi¬ 
tional kinematic variables in order to discriminate LQ signals from SM background processes and 
to enhance the signal-to-background ratio. The variables used are: 

- mu'. The dilepton invariant mass. 

- S't: The scalar sum of the transverse momenta of the two leading leptons and the two leading 
jets. 

- : The lowest reconstructed LQ mass in the event. The reconstructed masses of the two 

LQ candidates in the event and are defined as the invariant masses of the two 

lepton-jet pairs with the smallest difference (and < m™)- 

Signal regions are determined by optimising the statistical significance as defined in Ref. [56]. 
The optimisafion procedure is performed in a three-dimensional space constructed by mu, S t and 
for each of the signal mass points. Several adjacent mass points may be grouped into a 
single SR. The signal acceptance of the selection requirements is estimated to be w 50% in the 
fifijj channel and between 65% and 80% in the eejj channel (assuming yS = 1.0). The difference 
is due fo tighter quality selection requirements in the fifjjj channel used to prevent muon mis- 
measurements in MS regions with poor alignment or missing chambers. The optimised signal 
regions are presented in Table 1 together with the mass of the corresponding LQ hypothesis. Each 
LQ mass hypothesis is tested in only one signal region, where limits on cr x /? are extracted. 



LQ masses 
[GeV] 

niit 

[GeV] 

5t 

[GeV] 

[GeV] 

SRI 

300 

130 

460 

210 

SR2 

350 

160 

550 

250 

SR3 

400 

160 

590 

280 

SR4 

450 

160 

670 

370 

SR5 

500-550 

180 

760 

410 

SR6 

600-650 

180 

850 

490 

SR7 

700-750 

180 

950 

580 

SR8 

800-1300 

180 

1190 

610 


Table 1: The minimum values of , S t, and used to define each of the signal regions targeting different 
LQ masses in the eejj and pq/}' channels. Each signal region is valid for one or more mass hypotheses, as 
shown in the second column. 
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4.5 Background estimation 


The main SM background processes to the LQl and LQ2 searches are the production of Z/y*+jets 
events, tt events where both top quarks decay leptonically, and diboson events. Additional small 
contributions are expected from Z ^ tt and single-top processes. Multi-jets, VT-f-jets, tt (where 
one or more top quarks decays hadronically), and single-top events with mis-identified or non¬ 
prompt leptons arising from hadron decays or photon conversions can also contribute. These fake 
lepton backgrounds are estimated separately in the eejj and iijujj channels using the same data- 
driven techniques as described in Ref. [48] and are found to be negligible for the channel. 
Normalisation factors, derived using background-enriched control regions, are applied to the MC 
predictions for Z/y*-i-jets and ft backgrounds to predict as accurately as possible the background 
in the signal regions. These control regions are constructed to be mutually exclusive to the signal 
region and the assumption is made that normalisation factors and their associated uncertainties in 
the signal region are the same as in the background-enriched control regions. 


4.5.1 Control regions for Z/7*+jets and tt backgrounds 

Two control regions with negligible signal contributions are defined to validate the modelling ac¬ 
curacy of the MC simulated background events and to derive normalisation scale factors. The 
Z/y*-i-jets control region is defined by fhe pre-selecfion requiremenfs wifh an addifional require- 
menf of 60 < niee < 120 GeV in fhe eejj channel and 70 < <110 GeV in fhe ///ryy channel. 

These confrol regions define a pure sample of Z/y*-i-jefs evenfs. The tt confrol region is defined 
in bofh channels by applying fhe pre-selecfion requiremenfs, buf demanding exacfly one muon 
and one elecfron (bofh wifh pj above 40 GeV) in fhe offline selecfion insfead of fwo same-flavour 
lepfons. In fhe case of fhe tt confrol region for fhe eejj channel, fhe frigger requiremenf is modified 
by requiring a single isolafed elecfron wifh px above 24 GeV, which is fully efficienf relative fo fhe 
offline selecfion for elecfrons wifh px above 30 GeV. In bofh cases, fhe same selecfion criferia are 
applied fo dafa and MC evenfs. 

Normalisation factors are applied to fhe MC predicfions for fhe Z/y*-i-jels and tt background pro¬ 
cesses. They are obfained by performing a combined maximum likelihood til to fhe observed yields 
in fhe confrol regions and signal region under consideration. Sysfemafic uncerfainfies on fhe pre- 
dicfed MC yields relafed fo fhe uncerfainly on fhe cross-sections are faken info accounf by fhe fif 
fhrough fhe use of dedicated nuisance paramefers. The til procedure is performed using fhe HistFit- 
TER package [57], which is a fool based on fhe RooStats framework [58]. The normalisation scale 
factor obfained from a background-only fif for fhe Z/y*-i-jefs background in fhe eejj (ppjj) channel 
is 1.1 ± 0.2 (0.97 + 0.15), while fhe normalisation scale factor for ti is 1.10 + 0.05 (1.01 + 0.05). 
The filled background scale faclors have lillle sensilivily fo fhe inclusion of signal regions and fhe 
evenlual presence of a signal. 
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4.5.2 Kinematic distributions 


The distributions of the kinematic variables after performing the background-only fits in the control 
regions, and applying the event pre-selection requirements are shown in Figs. 1, 2 and 3 for the data, 
background estimates, and for three LQ masses of 300, 600 and 1000 GeV (with = 1.0). 



ITlee [GeV] 



[GeV] 


(a) eejj channel (b) mijj channel 

Figure 1; Distributions of the dilepton invariant mass (mu) in the eej] (left) and /i/i/}' (right) channels after 
applying the pre-selection cuts. The signal model assumes fi - 1.0. The last bin includes overflows. The 
ratio of the number of data events to the number of background events (and its statistical uncertainty) is also 
shown. The hashed bands represent all sources of statistical and systematic uncertainty on the background 
prediction. 


4.6 Systematic uncertainties 

The theoretical uncertainty on the NLO cross-section is taken into account for diboson, single¬ 
top, VT-i-jets, and Z ^ rr processes. For the two dominant backgrounds {tt and Z/y*-i-jets) the 
modelling uncertainties are estimated using the symmetrised deviation from unity of the ratio of 
data to MC events in the tt and Z/y*-i-jets control regions, which is fitted with a linear function for 
> 400 GeV. The modelling systematic uncertainty is then applied as a function of 5'T^ in the 
form of a weighting factor. The choice of 5 t for such a purpose is motivated by its sensitivity to 
mis-modelling of the kinematics of jets and leptons. It varies in the eejj (jiiijj) channel between 8 % 
(10%) and 25% (30%) for the Z/y*-i-jets background and between 6 % (10%) and 24% (40%) for 
the tt background. It increases for signal regions targeting higher /wlq. 

The jet energy scale (JES) uncertainty depends on pj and 77 and contains additional factors, which 
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(a) eejj channel 



(b) unij channel 


Figure 2; Distributions of the total scalar energy {S t) in the eejj (left) and (right) channels after applying 
the pre-selection cuts. The signal model assumes (3 - 1.0. The last bin includes overflows. The ratio of the 
number of data events to the number of background events (and its statistical uncertainty) is also shown. The 
hashed bands represent all sources of statistical and systematic uncertainty on the background prediction. 
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(a) eejj channel 


(b) mijj channel 


Figure 3: Distributions of the lowest reconstructed LQ mass ('«lq ) in the eejj (left) and p/p)' (right) channels 
after applying the pre-selection cuts. The signal model assumes (3- 1.0. The last bin includes overflows. The 
ratio of the number of data events to the number of background events (and its statistical uncertainty) is also 
shown. The hashed bands represent all sources of statistical and systematic uncertainty on the background 
prediction. 
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are used to correct for pile-up effects. They are derived as a function of the number of primary 
vertices in the event to take into account additional pp collisions in a recorded event (in-time pile- 
up), or as a function of the expected number of interactions per bunch-crossing to constrain past and 
future collisions affecting the measurement of energies in the current bunch-crossing (out-of-time 
pile-up). An additional uncertainty on the jet energy resolution (JER) is taken into account. The 
relative impact on the background event yields from the JES (JER) uncertainty is between 8% (1%) 
in SRI and 26% (1%) in SR8. The signal selection efficiency change due fo the JES uncertainties 
ranges between 3% in SRI and 1% in SR8, while the effect of the JER is negligible. 

The electron energy scale and resolution are corrected to provide better agreement between MC 
predictions and data. The uncertainties on these corrections are propagated through the analysis 
as sources of systematic uncertainty. Uncertainties are taken into account for the electron trig¬ 
ger (~ 0.1%), identification (~ 1%) and reconstruction (~ 1%) efficiencies, and for uncertainties 
associated with the isolation requirements (~ 0.1%). 

Scaling and smearing corrections are applied to the pj of the muons in order to minimise the differ¬ 
ences in resolution between data and MC simulated events. The uncertainty on these corrections is 
below 1%. Differences in fhe identification efficiency and in the efficiency of the trigger selection 
are taken into account and are less than 1%. 

QCD renormalisation and factorisation scales are varied by a factor of two to estimate the impact 
of higher orders on the signal production cross-section. The variation is found to be approximately 
14% for all mass points. The uncertainty on the signal cross-section related to the choice of PDE 
set is evaluated as the envelope of the prediction of 40 different CTEQ6.6 NEO error sets [24]. The 
uncertainty ranges from 18% at wilq = 300 GeV to 56% at wilq = 1300 GeV. These uncertainties are 
the same for all EQ generations. The effect on the choice of PDE set on the signal acceptance times 
reconstruction efficiency is estimated using the Hessian method [59]. The final PDE uncertainties 
on the signal samples are approximately 1% for most mass points, rising to 4% for some higher 
EQ masses. The impact of the choice of PDE set on the acceptance times reconstruction efficiency 
for each background process is estimated using the Hessian method (using the same method as for 
signals). The uncertainties range from 4% in the low-mass signal regions to 17% in the high-mass 
signal regions. 


4.7 Results 

The observed and expected yields in three representative signal regions for the eejj and the ppjj 
channels after the combined maximum likelihood fits are shown in Table 2 and Table 3, respect¬ 
ively. The fit maximizes the likelihood constructed using the two CRs and the SR under study. 
When contructing the likelihood, the signal stregth and the background scale factors are treated as 
free parameters, the systematic uncertainties are treated as nuisance parameters. 
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Yields 

1 eejj channel 

SRI 

SR6 

SR8 

Observed 

627 

8 

1 

Total SM 


11±2 

1.5 ±0.4 

Z/7* ^ U 

(3.2 ± 0.4) X 10" 

7±2 

1.3 ±0.4 

Z ^ TT 

2.1 ±0.3 

<0.01 

<0.01 

tt 

(2.4 ± 0.2) X 10" 

2.3 ± 0.5 

0.12 ±0.04 

Single top 

19 ±3 

<0.01 

<0.01 

Diboson 

22 ±3 

0.8 ±0.3 

<0.01 

Fake leptons 
(including IV+jets) 

34 ±6 

0.410 ±0.010 

0.033 ± 0.006 

mLQ = 300 GeV 


- 

- 

mLQ = 600 GeV 

nnmm 

231 ± 13 

- 

mLQ = 1000 GeV 


- 

5.2 ±0.3 


Table 2: Background and signal yields in three representative signal regions for LQs with masses mLQ = 300, 
600 and 1000 GeV for the eejj channel (assuming p - 1.0). The observed number of events is also shown. 
Statistical and systematic uncertainties are given. 


Yields 

1 fifijj channel 

SRI 

SR6 

SR8 

Observed 

426 

5 

1 

Total SM 


7.0 ± 1.2 

1.3 ±0.4 

Z/7* ^ a 

209 ± 18 

4.6 ± 1.0 

0.9 ± 0.3 

Z —> TT 

0.9 ±0.1 

< 0.01 

<0.01 

tt 

172 ± 18 

1.7 ±0.6 

0.18±0.11 

Single top 

14 ±5 

0.3 ± 0.4 

<0.01 

Diboson 

14 ±2 

0.5 ± 0.2 

0.19 ±0.05 

Fake leptons 
(including lV±jets) 

<0.01 

< 0.01 

<0.01 

mLQ = 300 GeV 


- 

- 

mLQ = 600 GeV 

■■■■ 

152 ± 18 

- 

mLQ = 1000 GeV 


- 

3.4 ± 1.3 


Table 3; Background and signal yields in three representative signal regions for LQs with masses mLQ = 300, 
600 and 1000 GeV for the ////// channel (assuming p - 1.0). The observed number of events is also shown. 
Statistical and systematic uncertainties are given. 
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No significant excess above the SM expectation is observed in any of the signal regions and a 
modified frequentist CL^ method [60] is used to set limits on the strength of the LQ signal, by 
constructing a profile likelihood ratio. Pseudo-experiments are used to determine the limits. 

The cross-section limits on scalar LQ pair-production are presented as a function of /? for both chan¬ 
nels in Fig. 4. Also shown are the results of the ATLAS searches for first- and second-generation 
LQs using 1.03 fb“^ data at ^/s= 1 TeV which also included searches in the evjj and //v/y de¬ 
cay channels and therefore provide better sensitivity at low values of /?. First (second)-generation 
scalar LQs are excluded for /3 = 1 at 95% CL for wilqi < 1050 GeV (/nLQ 2 < 1000 GeV). The 
expected exclusion ranges are the same as the observed ones. First (second)-generation scalar 
LQs are excluded for nrLQi < 650 GeV (mLQ 2 < 650 GeV) at /? = 0.2 and wilqi < 900 GeV 
('’^LQ 2 < 850 GeV) at yS = 0.5. 


5 Search for third-generation LQs in the bvrbv^ channel 

The ATLAS search for pair-production of third-generation supersymmetric partners of bottom 
quarks (sbottom, b) [18] is reinterpreted in terms of the LQ model, in the case where each LQ 
decays to a t7-quark and a Vr neutrino. In the original analysis, the b is assumed to decay via 
b^ bx^, and t via bx~ in the case where m^± - m^a is small and the x~ decay products are 
undetectable. The search is performed for final states with large missing transverse momentum 
with magnitude and two jets identified as originating from t^-quarks. The full ana¬ 

lysis strategy is covered in Ref. [18]. A complete description of the analysis, including treatment of 
systematic uncertainties on background processes can be found there, but the event selection and 
background estimation methods used are summarised here for clarity. 

5.1 Object and event selection 

Events are required to have exactly two fj-tagged [61] jets with pT > 20 GeV and \t]\ < 2.5, and E™*® 
> 150 GeV. Additional jets in the event are accepted if they have pT> 20 GeV and \r]\ <4.9. Events 
with one or more electrons (muons) with pi > 1 (6) GeV are vetoed. Candidate signal events 
are selected from data using a E™'®® trigger which is 99% efficient for events passing the offline 
selection. Several variables are defined and used to optimise the event selection: 

- A0min: The minimum azimuthal distance (A<p) between any of the leading three jets and the 

^miss 

Pj 

- niss'. The scalar sum of the pj of the leading two or three jets (depending on the signal 
region) and the 

- Ht: x The scalar sum of the px of all but the leading three jets. 


15 






mL02 [GeV] 


(d) 


Figure 4; The cross-section limits on scalar LQ pair-production times the square of the branching ratio as a 
function of mass (left) and the excluded branching ratio as a function of the LQ mass (right) to eq for the eej] 
channel (top) and to fj.q for the yUyuy channel (bottom). The +l(2)cr uncertainty bands on the expected limit 
represent all sources of systematic and statistical uncertainty. The expected NLO production cross-section 
(JS - 1.0) for scalar LQ pair-production and its corresponding theoretical uncertainty due to the choice of 
PDF set and renormalisation/factorisation scale are also included. The exclusion limits on LQl [12] and 
LQ2 [13] set by ATLAS in the eejj +evjj and fifijj +fivjj search channels using 1.03 fb“' of data collected at 
^[s -1 TeV are also shown. 
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- nihi,-. The invariant mass of the two Zj-tagged jets in the event. 

- mcT'- The contransverse mass [62], used to measure the masses of pair-produced heavy 
particles that decay semi-invisibly (i.e. decays where one of the decay products can be 
detected, but the other cannot). 

In the original analysis, different signal regions were optimised according to the masses of the 
third-generation squark and the lightest supersymmetric particle (LSP). In the case of the LQ 
model reinterpretation the signal regions corresponding to the case where the mass of the LSP 
is approximately zero have best sensitivity, but all the signal regions are retained for coherence 
with the original analysis. The different signal region definitions are given in Table 4. Signal re¬ 
gion A (SRA) has five differenf /tict fhresholds. Signal region B (SRB) is optimised towards fhe 
region where fhe squark and LSP masses are approximately equal. The signal region wifh fhe besf 
expecfed limif is used for each poinf in fhe exclusion plofs. 


5.2 Background estimation 

The dominanf background process is fhe production of Z bosons in associafion wifh heavy-flavour 
Jefs where fhe Z boson subsequenfly decays to fwo neufrinos (Z(—> yy) -i- bb). Ifs confribufion is 
esfimafed from dafa in an opposife-sign dilepfon confrol region. Top quark pair-production {tt) 
and W bosons produced in association wifh heavy flavour quarks also confribufe significanfly and 
are normalised in dedicated confrol regions before being exfrapolaled fo fhe signal regions using 
MC simulation. Differenf confrol regions are defined for each signal region, requiring one or fwo 
leptons plus addifional requiremenfs similar fo fhe corresponding signal region. The confribufions 
from Z-i-jels, IT-f-jels, and lop quark production are esfimafed simullaneously wifh a profile likeli¬ 
hood fil to fhe fhree confrol regions. Confribufions from diboson and tt+WjZ processes are estim¬ 
ated from MC simulalion in all regions. The confribufion from mulfi-jel evenls is esfimafed from 
dafa by faking well-measured mulfi-jel evenls from dafa and smearing fhe jefs wifh jel response 
functions laken from MC simulalion and validaled in dafa. This procedure is described in delail in 
Ref. [63]. The confribufion from mulfi-jel evenls in signal regions is found fo be negligible. 

5.3 Results 

The number of dafa evenls observed in each signal region is reported in Table 5, logelher wifh 
fhe SM background expecfafion affer fhe background-only fif, and fhe expecfed number of signal 
evenls for differenf LQ masses. The signal accepfance efficiency is around 2% for all buf fhe lowesl 
LQ masses largeled (dropping to 0.27% efficiency for wilq = 200 GeV). All sources of syslemalic 
and slalislical uncertainly are laken info accounl. The dominanf systematic uncertainlies on fhe 
background prediclion are fhe jel energy scale (JES: 1-5%) and resolufion (JER: 1-8%), and fhe 
fj-fagging uncerfainfy (2-10%). Defeclor-relafed systematic uncerlainlies on fhe signal prediction 
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Description 

Signal Regions 

SRA 

SRB 

Event cleaning 

Common to all SR 

Lepton veto 

No e/p after overlap removal with px > 7(6) GeV for e(p) 

^miss 

> 150GeV 

> 250 GeV 

Leading jet pj 

> 130GeV 

> 150 GeV 

Second jet px 

>50 GeV 

> 30 GeV 

Third jet px 

veto if > 50 GeV 

> 30 GeV 

lead jet) 

- 

> 2.5 

Zi-tagging 

leading 2 jets 

(px>50GeV, |p|<2.5) 

2nd- and 3rd-leading jets 

(px > 30GeV,|77| <2.5) 


A^rain 

> 0.4 

>0.4 

(OTeff = £'t'“ + Pt + Pj) 

£7“/meff > 0.25 

- 

(meff = + p!^ + p" + p") 

- 

/meff > 0.25 

mcT 

> 150, 200, 250, 300, 350 GeV 

- 

Ht.3 

- 

< 50 GeV 

nibb 

> 200 GeV 

- 


Table 4: Summary of the event selection in each signal region for the bvrbvT channel [18]. 
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SRA, nicT > 



SRB 


150 GeV 

200 GeV 

250 GeV 

300 GeV 

350 GeV 


Observed 

102 

48 

14 

7 

3 

65 

Total SM 

94± 13 

39 ±6 

16 ±3 

5.9 ± 1.1 

2.5 ± 0.6 

64± 10 

Top quark 

11.1 ± 1.8 

2.4 ± 1.4 

0.4 ± 0.3 

<0.01 

<0.01 

41 ±7 

Z production 

66± 11 

28 ±5 

11 ±2 

4.7 ± 0.9 

1.9 ±0.4 

13 ±4 

W production 

13±6 

5±3 

2.1 ± 1.1 

1.0 ±0.5 

0.5 ±0.3 

8±5 

Others 

4.3 ± 1.5 

3.4 ± 1.3 

1.8 ±0.6 

0.12±0.11 

0 

'^•^^-0.10 

2.0 ± 1.0 

Multi-jet 

0.2 ± 0.2 

0.06 ± 0.06 

0.02 ± 0.02 

< 0.01 

< 0.01 

0.16 ±0.16 

otlq = 300 GeV 

(8.5 ± 0.2) X 10^ 

435 ± 17 

96 ±8 

7±2 

0.6 ± 0.6 

68 ±7 

“LQ = 600 GeV 

21.9 ±0.4 

19.0 ±0.4 

15.6 ±0.4 

12.0 ±0.3 

8.7 ±0.3 

1.8±0.1 


Table 5: For each signal region in the bv-jjVr channel, the observed event yield is compared with the back¬ 
ground prediction obtained from the fit. Signal yields for different values of tolq (assuming/? = 0.0) are given 
for comparison. The category ‘Others’ includes the diboson and tt+WIZ processes. Statistical, detector- 
related and theoretical systematic uncertainties are included, taking into account correlations [18]. 


are dominated by uncertainties on the Z^-tagging efficiency (~30%). The second-largest source of 
uncertainty is due to the JES and is around 3%. 

The uncertainties on the signal production cross-section are estimated using the methods described 
in Section 4.6. These uncertainties are the same for all LQ genertions but the uncertainty due to the 
choice of PDF set varies with wtlq. Since the third-generation analyses consider a different mass 
range to the first- and second-generation analyses, in this case the uncertainty due to the choice of 
PDF set ranges from 7.1% at wlq = 200 GeV to 30% at mLQ = 800 GeV. Effects on the acceptance 
due to the choice of PDF set are negligible. 

No significant excess above the SM expectation is observed in any of the signal regions. Figure 5 
shows the observed and expected exclusion limits for the scalar FQ3 pair-production scenario ob¬ 
tained by taking, for each signal mass configuration, the signal region with the best expected limit. 
These limits are obtained using the methods described in Section 4.7. These methods compare 
the observed numbers of events in the signal regions with the fitted background expectation and 
accounting for signal contamination in the corresponding CRs for a given model. Pair-produced 
third-generation scalar FQs decaying to bvrbVr are excluded at 95% CF for /«lq 3 < 625 GeV. The 
expected excluded range is mLQ 3 < 640 GeV. 


6 Search for third-generation LQs in the tVrtVr channel 

The ATFAS search for pair-production of the supersymmetric partner of the top quark (stop quark, 
t) [19] is reinterpreted in terms of the FQ model, in the case where each FQ decays to a top quark 
and a Vr neutrino. The original analysis has dedicated signal regions targeting t decays into tx^ and 
the subsequent semileptonic decay of the ti pair. Events compatible with ti plus extra are se¬ 
lected with final states containing one isolated lepton, jets, and A complete description of the 
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mLQ3 [GeV] 


Figure 5: The expected (dashed) and observed (solid) 95% CL upper limits on third-generation scalar LQ 
pair-production cross-section times the square of the branching ratio to bvj as a function of LQ mass, for the 
bvrbvj channel. The +l(2)cr uncertainty bands on the expected limit represent all sources of systematic and 
statistical uncertainty. The expected NLO production cross-section (JS - 0.0) for scalar LQ pair-production 
and its corresponding theoretical uncertainty due to the choice of PDF set and renormalisation/factorisation 
scale are also included. 
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analysis strategy, including the treatment of systematic uncertainties on background processes can 
be found in Ref. [19]. The event selection and the background estimation methods are summarised 
here for clarity. 

6.1 Object and event selection 

Events are required to contain exactly one electron with pj > 25 GeV and \r]\ < 2.47, or muon 
with pt > 25 GeV and \r]\ < 2.4. Events containing more than one electron or muon with looser 
identification and pT requirements (10 GeV for both) are vetoed. In some signal regions, events are 
vetoed if they are consistent with containing a hadronically decaying t lepton. Events are required 
to have a minimum of four jets with pj > 20 GeV and \r]\ < 2.5, with at least one of these passing 
fj-tagging requirements [61]. In addition, selected events must have > 100GeV. Several 
variables are used to further select signal events and reject background processes: 

- mj\ The transverse mass of the electron or muon and the E™®®. 

- aniT and These are two variants on the stransverse mass (/nT, 2 ) [64-66] which is a 

generalisation of the transverse mass when applied to signatures with two invisible particles 
in the final sfafe. The asymmefric sfransverse mass amx, aims fo rejecf dilepfonic tt evenfs 
where one of fhe lepfons is nol reconsfrucfed or is oufside fhe accepfance (and Iherefore adds 
fo fhe E™*® of fhe even!). The second implemenfafion of Ibis variable, fhe t sfransverse mass 

fargefs tt events where one top decays leptonically and the other top decays into a r 
that subsequently decays hadronically. 

- topness: This variable is designed to reject dileptonic tt events where one lepton is assumed 
to be lost, as detailed in Ref. [67]. The topness variable is based on the minimisation of a 
;t'^-type function. 

- mhad-top: This quantity is used to reject dileptonic tt events but retain signal events that con¬ 
tain a hadronically decaying on-shell top quark, as in the EQ t+Vr and t^ —> tk? scenarios. 

- A0(jet[ The azimuthal opening angle between the leading or sub-leading jet and 

used to suppress multi-jet events where is aligned with one of the leading two 

jets. 

- ! V^‘ An approximation of the E™*® significance, where T/j is defined as fhe scalar 
Pt sum of fhe leading four jefs. 

- An objecf-based missing transverse momentum, divided by the per-event resolution 
of the jets, and shifted to the scale of the background [68]. 
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The variables listed above are used to define three cut-and-count SRs and one shape-fit SR. Table 6 de¬ 
tails the event selections for these signal regions. The SR labelled tN_boost targets LQ/stop 
masses of >700 GeV and takes advantage of the ‘boosted’ topology of such a heavy parent particle. 
The selection assumes that either all decay products of the hadronically decaying top quark, or at 
least the decay products of the hadronically decaying W boson, collimate into a jet reconstructed 
with a radius parameter R = 1.0 [69,70]. 


6.2 Background estimation 

The dominant sources of background are the production of tt events and VT-i-jets where the W 
boson decays leptonically. Other background processes considered are single top, dibosons, Z-i-jets, 
tt produced in association with a vector boson (ttV), and multi-jets. 

The predicted numbers of ti and IT-i-jets background events in the SRs are estimated from data using 
a fit to the number of observed events in dedicated control regions. Each SR has an associated CR 
for each of the tt and IT-i-jets backgrounds. The CRs are designed to select events as similar as 
possible to those selected by the corresponding SR while keeping the contamination from other 
backgrounds and potential signal low. This is achieved by e.g. requiring that 60 < ni'Y < 90 
GeV and in the case of the fT-i-jets CR, inverting the ft-jet requirement so that it becomes a Z^-jet 
veto. The simulation is used to extrapolate the background predictions into the signal region. The 
background fit predictions are validated using dedicated event samples, referred to as validation 
regions (VRs), and one or more VR is defined for each of these. Most VRs are defined by changing 
the mx windows to 90 < mj <120 GeV. The VRs are designed to be kinematically close to the 
associated SRs to test the background estimates in regions of phase space as similar as possible to 
the SRs. 

The multi-jet background is estimated from data using a matrix method described in Refs. [71,72]. 
The contribution is found to be negligible. All other (small) backgrounds are determined entirely 
from simulation and normalised to the most accurate theoretical cross-sections available. 

6.3 Results 

The number of events observed in each signal region is reported in Tables 7 and 8, together with the 
SM background expectation and the expected number of signal events for different LQ masses. The 
signal acceptance is between 1.5% and 3% depending on the LQ mass. All sources of systematic 
uncertainty and statistical uncertainty are taken into account. The dominant sources of uncertainty 
on the background prediction come from uncertainties related to the JES, JER, ttbackground mod¬ 
elling, the fi-tagging efficiency, and statistical uncertainties. 
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tN_diag 

tN_med 

tN_high 

tN_boost 

Lepton 

= 1 lepton 

Jets 

> 4 with pt > 

60,60,40,25 GeV 

> 4 with pt > 

80,60,40,25 GeV 

> 4 with Pt > 

100, 80,40,25 GeV 

> 4 with Pt > 

75,65,40, 25 GeV 

i’-tagging 

> Ifc-tag (70% elf.) amongst four selected jets 

Large-R jet 

- 

> l,pT> 270GeV 

and m > 75 GeV 

A0(jet!,“®‘='',/J“‘“) 

- 

>0.85 

^miss 

> 100GeV 

> 200 GeV 

> 320 GeV 

> 315 GeV 

nij 

> 60 GeV 

> 140 GeV 

> 200 GeV 

> 175 GeV 

amj 

- 

> 170 GeV 

> 170 GeV 

> 145 GeV 


- 

- 

> 120 GeV 

- 

Topness 

- 

- 

- 

> 7 

^had-top 

6 [130, 205] GeV 

6 [130, 195] GeV 

6 [130, 250] GeV 


r-veto 

loose T particle ID 

- 

- 

modified, see [19]. 

Ai?(Z7-jet, f) 

< 2.5 

- 

< 3 

< 2.6 


> 5 GeV'/2 

- 

r/miss 

llrj. ■ 

T,sig 

- 

> 12.5 

> 10 


> 0.8 (i = 1,2) 

> 0.8 (/ = 2) 

- 

> 0.5,0.3 (( = 1,2) 

Exclusion setup 

shape-fit in ott and 

^miss 

cut-and-count 


Table 6: Selection criteria for the four SRs (tN_diag, tN_med, tN_high, and tN_boost) employed to search 
for LQ3 LQ3 ^ fv^rtv events [19]. The details of the limit-setting procedure for the exclusion setup can be 
found in Section 6.3. 


23 




































tN med 

tN high 

tN boost 

Observed 

12 

5 

5 

Total SM 

13 ± 2 

5.0 ± 0.9 

3.3 ± 0.7 

tt 

6.5 ± 1.7 

2.0 ± 0.6 

1.1 ±0.4 

W+jets 

2.1 ±0.5 

0.9 ± 0.3 

0.28 ±0.14 

Single top 

1.1 ±0.5 

0.54 ±0.19 

0.39 ±0.15 

Diboson 

1.4 ± 0.6 

0.9 ± 0.3 

0.7 ± 0.3 

Z+jets 

0.009 ± 0.005 

0.003 ± 0.002 

0.004 ± 0.002 

ttV 

2.0 ± 0.6 

0.8 ± 0.3 

0.9 ± 0.3 

mLQ = 300 GeV 

20 ±3 

3.4 ± 1.1 

3.8 ± 1.2 

niLQ = 600 GeV 

10.7 ± 0.3 

7.9 ± 0.3 

8.9 ± 0.3 


Table 7: The number of observed events in the three cut-and-count signal regions, together with the expected 
number of background events and signal events for different LQ masses (assuming fi - 0.0) in the tVrivr 
channel [19]. 


Detector-related systematic effects are evaluated for signal using the same methods used for the 
backgrounds (see Ref. [19] for details). The dominant detector-related systematic effects are the 
uncertainties on the JES (4%) and the Z^-tagging efficiency (3%). 

The uncertainties on the signal production cross-section are estimated using the methods described 
in Section 4.6. The effect on the choice of PDF set on the signal acceptance is less than 1% for 
most mass points, but increases to 1.7% for ^lq = 800 GeV. 

Similar methods as described in Section 4.7 are used to assess the compatibility of the SM background- 
only hypothesis with the observations in the signal regions. The observed number of events is found 
to agree well with the expected number of background events in all signal regions. No significant 
excess over the expected background from SM processes is observed and the data are used to derive 
one-sided limits at 95% CL. The results are obtained from a profile likelihood-ratio test following 
the CL^ prescription [60]. The likelihood of the simultaneous fit is configured to include all CRs 
and one SR or shape-fit bin. The ‘exclusion setup’ event selection is applied (see Table 6), and all 
uncertainties except the theoretical signal uncertainty are included in the fit. 

Exclusion limits are obtained by selecting a priori the signal region with the lowest expected CL^ 
value for each signal grid point. The expected and observed limits on the LQ3 LQ3 —> tVrrtv pro¬ 
cess are shown in Fig. 6. Third-generation scalar LQs decaying to tv-^tYr are excluded at 95% CL 
in the mass range 210 < mLQ 3 < 640 GeV. The expected exclusion range is 200 < mLQ 3 < 685 GeV. 
The limits for stop production in the case where the neutralino is massless are slightly stronger than 
the limits set on LQ3 production since the nominal stop limits consider a mostly right-handed stop. 
This leads to the top quarks being polarised in such a way that the acceptance increases. The limit 
worsens at low mass, due to the effect of greater contamination from top backgrounds. 
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Figure 6: The expected (dashed) and observed (solid) 95% CL upper limits on the third-generation scalar 
LQ pair-production cross-section times the square of the branching ratio to fVr as a function of LQ mass, 
for the fVrtVr channel. The ±\{2)cr uncertainty bands on the expected limit represent all sources of system¬ 
atic and statistical uncertainty. The expected NLO production cross-section {[i - 0.0) for third-generation 
scalar LQ pair-production and its corresponding theoretical uncertainty due to the choice of PDF set and 
renormalisation/factorisation scale are also included. 
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125 < < 150 GeV 

120 < niT < 140 GeV 

tN_diag 
125 < L™'*" < 150 GeV 
niT > 140 GeV 

> 150 GeV 
120 < niT < 140 GeV 

> 150 GeV 
Tot > 140 GeV 

Observed 

117 

163 

101 

217 

Total SM 

(1.4 ± 0.2) X 10^ 

(1.5 ± 0.2) X 10^ 

98 ± 13 

(2.4 ± 0.3) X 10^ 

tt 

(1.2 ± 0.2) X 10- 

(1.4 ± 0.2) X 10^ 

85 ± 12 

(2.0 ± 0.3) X 10^ 

Vf-l-jsts 

7±3 

6±3 

4.6 ± 1.5 

10 ±4 

Single top 

5±2 

6± 2 

6± 2 

9±4 

Diboson 

0.29 ±0.18 

0.8 ± 0.5 

0.3 ± 0.3 

0.30 ±0.15 

Z-l-jets 

0.17 ±0.08 

0.24 ±0.12 

0.30 ±0.15 

0.5 ± 0.3 

ttV 

1.5 ±0.5 

2.9 ± 0.9 

2.5 ± 0.8 

11 ±3 

mLQ = 300 GeV 

28 ±3 

77 ± 6 

64 ±5 

269 ± 10 

niLQ = 600 GeV 

0.15 ±0.04 

0.62 ± 0.08 

0.83 ± 0.09 

18.8 ± 0.4 


Table 8; The number of observed events in the shape-fit signal region, together with the expected number of 
background events and signal events for different LQ masses (assuming p — 0.0) in the fv^fvv channel [19]. 


7 Summary and conclusions 


Searches for pair-production of first-, second- and third-generation scalar leptoquarks have been 
performed with the ATLAS detector at the LHC using an integrated luminosity of 20 fb“^ of data 
from pp collisions at = 8 TeV. No significant excess above the SM background expectation is 
observed in any channel. The results are summarised in Table 9. 


Decay channel 

Excluded ran 
Expected 

ge (95% CL) 

Observed 

eejj{fi= 1.0) 

mLQi < 1050 GeV 

mLQi < 1050 GeV 

mj(J3= TO) 

inLQi < 1000 GeV 

™lq 2 < 1000 GeV 

bvrbvr(J3 = 0.0) 

tnLQ 3 < 640 GeV 

"!lq 3 < 625 GeV 

tvpvrip = 0.0) 

200 < mLQ 3 < 685 GeV 

210 < OTLQ 3 < 640 GeV 


Table 9; Expected and observed exclusion ranges at 95% CL for each of the four LQ decay channels con¬ 
sidered. 


The results presented here significantly extend the sensitivity in mass compared to previous searches. 
Low-mass regions are also considered and limits on the cross-sections are provided for the different 
final sfafes analysed. Since /? is nol consfrained by fhe fheory, searches in fhe low mass regions are 
also imporfanf in order fo exfracf limifs for low-/? values for fhe LQl and LQ2 analyses. 
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